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Abstract 



We give an exact spherically symmetric solution for the Einstein-scalar 
field system. The solution may be interpreted as an inhomogeneous dynamical 
scalar field cosmology. The spacetime has a timelike conformal Killing vector 
field and is asymptotically conformally flat. It also has black or white hole- 
like regions containing trapped surfaces. We describe the properties of the 
apparent horizon and comment on the relevance of the solution to the recently 
discovered critical behaviour in scalar field collapse. 
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The gravitational collapse of distributions of matter is one of the most important research 
areas in general relativity. The essential question posed is whether and under what initial 
conditions a black hole or naked singularity forms in the collapse. One of the aims of studying 
such problems is to test the cosmic censorship conjecture one form of which states that 
gravitational collapse produces black holes. 

A model system for studying this question is provided by the Einstein equations mini- 
mally coupled to a massless scalar field. While this full system appears to be intractable, 
the simplified set of equations obtained by imposing spherical symmetry is easier to han- 
dle. Without any matter fields, the spherically symmetric metric does not contain any field 
degrees of freedom. Therefore, with a scalar field, the system is effectively a two dimen- 
sional field theory and it can be described by a single two dimensional nonlinear differential 
equation 

There are a number of exact solutions known for this system, almost all of which are 
either static or depend only upon the time coordinate The first non-static solutions has 
been given by Roberts |]], (which are different from the one we give below). The equations 
have been studied in detail by Christodoulou 0] who established, among other things, that 
there exist regular solutions for arbitrarily long times for particular types of initial data. 

The model has also been studied numerically and there are a number of interesting 
numerical results. The first results obtained by Goldwirth and Piran indicated that 
there is a class of initial data that leads to black hole formation. More recently it has 
been shown by Choptuik || that, for large classes of initial data, there is critical behaviour 
at the onset of black hole formation: the black hole mass Mbh is given by the equation 
Mbh — K\c — c*| 7 , where K is a constant, c is any one of the parameters in the initial 
data for the scalar field, c* is a critical value of the parameter, and 7 ~ .37 is a universal 
exponent. The remarkable feature of this result is that 7 appears to be independent of a 
set of particular shapes of the initial data, and is universal in this sense. It has been shown 
by Abraham and Evans || that the same mass equation is obtained for the axisymmetric 
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collapse of gravitational radiation. Thus the critical behaviour appears to be independent 
of not only the type of matter fields, but also the symmetries of the system. 

It would be very useful to understand the universality of this result analytically. A 
modest approach is to attempt to find an exact solution describing scalar field collapse and 
to see if one can read off the critical behaviour by calculating the mass of the black hole. 

Here we describe an exact solution for scalar field collapse and discuss some of its prop- 
erties. While the solution we present does not describe a realistic collapse corresponding 
to the classes of initial data used in the numerical work mentioned above, it appears to be 
among the few exact non-static solutions known for this system. 

The Einstein-scalar field equations we consider (in units G — c — 1 ) are 





(1) 



which may be written in the form 



(2) 



(where the subscript on denotes partial differentiation). 



The spherically symmetric solution we obtain is 



ds 2 = (at + b) {-f{r)dt 2 + r\r)dr 2 ) + R 2 (r,t)(d6 2 + sin 2 9dcj> 2 ), 



(3) 



where 



fir) 



(1 - 2c/r) 



a 



R 2 (r-t) 



(at + b)r 2 (l-2c/r) 



(4) 



The scalar field is 





where a, b, c, d are constants, 



a = ±\/3/2, and the overall sign of is independent 
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of a. The coordinate ranges and the values of the constants are coupled (for the metric 
to be Lorentzian): — b/a <t<oo and 2c < r < oo (c > 0). 

We note that for a ^ there is a coordinate transformation, that eliminates the constant 
b. However, since we will also be interested in the metrics for which a = 0, it is useful to 
keep the general form (^). We note also that when 6 = 0, there are two ranges for t, namely 
< t < oo for a > 0, and — oo < t < for a < 0. As shown below these correspond to white 
and black hole-like solutions respectively, and simply reflect the choice of the arrow of time. 

When a ^ 0, the coordinate transformation T = at + b does not eliminate 

a in the metric, which is an ovserall scale. The constant d in the 

scalar field is a trivial additive constant. 

The only Killing vectors of the metric (3) are the three associated with spherical sym- 
metry. The metric also has the conformal Killing vector field V = d/dt such that 

CyQr = ^9^. (6) 

Asymptotically (r — > oo), the metric is conformal to the Minkowski metric. The locus of 
points at r = 2 is a timelike curvature 

singularity, and so the 'horizon' is shrunk to a point as in the static metric given by 
Janis, Newman and Winicour (JNW), and others ||||. There is also a spacelike singularity 
at t = —b/a. The solution may be interpreted as a scalar field cosmology, since it is not 
asymptotically flat. 

For the JNW metric, the functions f(r) have arbitrary exponents specified by integration 
constants rather than the fixed a = ±v3/2 above. Thus our metric is conformal to the JNW 
metric with this exponent fixed, and with (at + b) as the conformal factor. In fact for a = 
we recover one of the JNW metrics. Therefore a distinguishes static from 

non-static metrics. This is similar to the parameter k in 

Friedman-Robertson- Walker cosmologies, which discretely distinguishes 

the spatial curvatures of the metrics. 

The case c = gives a homogeneous time dependent solution. This is equivalent to the 
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r — > oo limit of (|3|). The parameter c therefore distinguishes homogeneous from inhomo- 

geneus time dependent solutions. 

For comparison with recent numerical work |||9|], where the scalar 
field and its time derivative are specified as part of the initial data on 
a spacelike hypersurface, we note that this data for our solution is 

0(r, t = t ) = — ^= ln(d (1-2 c/r) v% ( a t + d) ±V *), 

The asymptotic behaviour of is 

0( r = oo,t = t ) = -^ln(d(at + b) ±VTi ). (8) 
4 V 71 " 

These are not the initial data associated with the standard collapse situation ||||, where 
the data is typically an ingoing pulse 

with a specified amplitude and width. 

The Ricci scalar is 

= \2ca\r - c) - 3a 2 r 2 _ 2_c 2c 2 (l - a 2 ) _ 2c 2+a 

2r 2 {at + bf 1 r ] (at + 6)r 4 1 r 1 K 1 

which shows that curvature singularities are present at r = 2c, and at t — —b/a. 

Since the metric is not static, it is of interest to investigate the existence and properties 
of the apparent horizon and to calculate the mass function. The apparent horizon is the 
3-surface on which outgoing or ingoing null rays are momentarily stationary. The presence 
of the horizon for asymptotically flat spacetimes indicates that there is a black or white hole, 
and a measure of the mass within it is given by the mass function evaluated at the horizon. 
This measure of the black hole mass has been used in recent numerical work on the collapse 
of ingoing matter pulses f|,(|. Although the metric (0) is not asymptotically flat, we can 
nevertheless see whether there are horizons. 

The apparent horizon surface is given by 
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g aP R, a R,p=0, (10) 



which for our metric gives the equation 

a 2 



-[r - c(l + «)](!- 2 c/r)"" 1 (11) 



at ah + b r 2 

This equation has no non-trivial solution for a = which corresponds to the static JNW 
metric, wheras for a ^ there is always an evolving apparent horizon. 

The apparent horizon can in general be spacelike, null or timelike in different spacetime 
regions. This is easily determined by calculating the ratio of the slopes of the apparent 
horizon and the outgoing null ray. This ratio for our metric (with a ^ 0) is 

t AH,r_ l (l-2c/r) 

t N>r 2[l-c(l + a)/rf 1 ' 

The second term on the right hand side is always positive, therefore the apparent horizon is 
spacelike for all r > 2c. It is null only at r = 2c. 

The scalar field is not singular at the apparent horizon. 

Figures 1 and 2 are plots of the horizon for a = ±y/3/2 and c = 1,6 = 0. The main 
features remain unaltered for all values of 6, c. The horizon forms at r = 2c, t — and 
grows in size forever. For a = y/3/2, the light cones are collapsed to a vertical line at r = 2c 
and open up to a slope of ±l/6 2 as r — > oo, whereas for a = — Vo/2 the cones collapse to a 
horizontal line at r = 2c. 

It is of interest to note a number of other features of the apparent horizon. By computing 
the expansions of the spacelike symmetry 2-spheres 

ds 2 = R 2 {r,t){d6 2 + sm 2 6d<j) 2 ) (13) 

along future pointing null directions orthogonal to the spheres, one can determine whether 
the apparent horizon is past or future, and inner or outer, and what region is trapped (see 
pTjj] for a general discussion of apparent horizons). 

The expansions 8± of the area 2- form u = R 2 (r, t)sin# d6 A dcf) of the 2-spheres are 
defined by 



£ 1± lu = ± w, (14) 
where £ denotes the Lie derivative and 

* + :=£ + ff i--^-/ 2 ^ (15) 
ot Or ot Or 

are the outgoing and ingoing future pointing null directions. With 6 = and a > 0, (so that 
< t < oo), we find 

e± = 7 ± t 1 -, (i6) 

with as given in Eq. (11). Thus it is the ingoing expansion 6- that vanishes at the 
apparent horizon, while the outgoing expansion is 6+ = 2/tAH > 0. This implies that the 
horizon is a past horizon. For a given value of r, the symmetry 2-spheres are trapped surfaces 
for t < tjyjj, since this is the region where both the ingoing and outgoing light rays have 
positive expansions. Similarly, the region t > tAH is a normal region where the outgoing 
light expansion (#+) is positive and the ingoing one (8 '_) is negative. 

We note also that if £i + 6_\ah < the horizon is an outer one (otherwise it is inner). 



For the metric (Q), we find 

l AH l N,r 

From flTJD it follows that this is always less than zero except at the singularity where it is 
zero. 

Summarizing the above results, the apparent horizon is a past outer one and is spatial 
everywhere except at r = 2c, where it is null. The scalar field flows from the past trapped 
region t < tAH (white hole) into the untrapped region t > tAH- An observer in the untrapped 
region sees both the t — initial singularity and the one at r = 2c. 

When 6 = 0, there is a change for the time reversed case which corresponds to a < 
and — oo < t < 0. The horizon is now a future outer one and corresponds to a black hole 
situation. The future singularity at t = is covered by the spacelike horizon and the region 
is a black hole. 



We now turn to a discussion of the mass function defined by 

m(r,t) = ^(l-g^R, a R, (3 ). (18) 

This function for static or stationary asymptotically flat spacetimes gives the ADM mass in 
the asymptotic limit. Its value on the apparent horizon is obtained by substituting tAH(f) 



from (|TT[) into 



Rah , \ i \a\r 2 {l-2c/r) l - a 
M , = ^L = m (r,t)\ AH = \ V f ' ' • (19) 

2 V 8 Jr-c(l + a) 

The mass M is zero at r = 2c and grows as r 3//2 for large r. We note that it is independent 
of b. For a = 0, M = 0, which corresponds to the static JNW solution. This corresponds to 



the fact that there is no solution to the apparent horizon equation (|TT|) for a = 0. 

Superficially, the mass function may be viewed as describing a form of critical behaviour 
since for fixed r we can write M = constant(a — a*) 1 / 2 , with the 'critical' value a* = 0, 
and with one half as the 'exponent'. Similarly, one can do this for the parameter c in (|T9|). 
which gives a different exponent. This, however, does not shed light on the collapse of initial 
pulses of scalar field, since the initial data (0) does not correspond to this situation. 

The metric is nevertheless an exact solution and may be viewed as a model in which the 
formation and evolution of an apparent horizon can be studied exactly. The value of the 
parameter a in the solution determines whether an apparent horizon forms. It would be of 
much interest to find an exact solution corresponding to a realistic collapse, as this would 
provide an analytical model for viewing the scaling and critical behaviour obtained in recent 
numerical work [§,0. 
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FIGURE CAPTIONS 

Figure 1: The apparent horizon for a = anc ^ c = 1,6 = 0. The light cones 

illustrate the spacelike character of the horizon. The singularities are at 
r = 2c and t — 0. The trapped region is t < tAH- 

Figure 2: The apparent horizon for a = — ^ and c = 1, b = 0. 
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